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Phase diagram of hot and dense QCD constrained by the Statistical Model 
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We propose a prescription to constrain the chiral effective model approach to the QCD phase 
diagram using the thermal Statistical Model, which is a description consistent with the experimental 
data at the freeze-out. In the transition region where thermal quantities of hadrons blow up, 
deconfined quarks and gluons should smoothly take over the degrees of freedom from hadrons in 
the Statistical Model. We use the Polyakov-loop coupled Nambu-Jona-Lasinio (PNJL) model as 
an effective description in the quark side. We require that the validity regions of these descriptions 
should have an overlap on the phase diagram, which constrains model uncertainty. Our results 
favor a phase diagram with the chiral phase transition located at slightly higher temperature than 
deconfinement . 
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Introduction Exploration of the QCD (Quantum 
Chromodynamics) phase diagram, particularly toward 
higher baryon-density regime, is of increasing importance 
in both theoretical and experimental sides ^1]. From the 
theoretical point of view, so far, only the lattice-QCD 
simulation jl], Q is the first-principle approach at work 
to the QCD phase transitions — chiral restoration and 
quark deconfinement. The functional renormalization 
group method is also developing as a promising non- 
perturbative tool [3j. The chiral condensate (ij^ip) and 
the Polyakov loop $ are the (approximate) order param- 
eters for chiral restoration and quark deconfinement, re- 
spectively, which are gauge invariant and measurable on 
the lattice. The lattice-QCD simulation is, however, of 
practical use only when the baryon chemical potential 
fiB is sufficiently smaller than the temperature T. For 
Hb/T > 1 the notorious sign problem prevents us from 
extracting any reliable information from the lattice-QCD 
data [1, 4]. 

The effective model study is an alternative and prag- 
matic approach toward the phase diagram of dense QCD. 
The idea is the following; one starts with some models 
that yield a reasonable description of hadron properties 
in the vacuum and then puts them in a finite-T and/or 
finite- /iB environment. What is recognized nowadays as 
the "QCD phase diagram" is actually a theoretical con- 
jecture based on various effective model studies. 

Along this line the Polyakov-loop coupled chiral 
models such as the PNJL (Polyakov-Nambu-Jona- 
Lasinio) 0,11] and the PQM (Polyakov-Quark- Meson) 0, 
Q models are successful to handle (ipip) and $ on the 
equal footing. Besides, the Polyakov loop potential 
is determined by $ and the pressure p measured in the 
lattice simulation of the pure gluonic theory. This means 
that the model includes the pressure contribution from 
gluons as well as quarks, so that the model is able to 
deal with the full thermodynamics comparable with the 
full lattice-QCD simulation. The point is that the dy- 
namics of transverse gluons A[ is under the control of 
the deconfinement order parameter $ and thus is to be 



encompassed in the Polyakov loop potential U[^], while 
the Polyakov loop itself is expressed in terms of the lon- 
gitudinal gluon Ai. 

Since theory instruments to examine both (ipip) and $ 
are now in our hands, it is intriguing to address the fol- 
lowing question; whether the chiral and deconfinement 
phase transitions would go on simultaneously or sepa- 
rate after all when the baryon density increases. There 
are then two key issues. One is the so-called QCD (chi- 
ral) critical point (which is often called the critical end- 
point) at which the chiral and the baryon number sus- 
ceptibilities diverge JoMlII and the higher moments are 
even more singular [12| . The other one is a triple-point- 
like region associated with the appearance of quarkyonic 
matter 13, [l6| where the baryon abundance surpasses 
mesons. 

One reasonable way to characterize quarkyonic matter 
for finite- iVc QCD is to use two order parameters $ = 
and the quark (baryon) number density (tp'^ip) ^ 0, which 
would definitely work for = oo |13| . In principle this 
statement is not directly related to chiral symmetry, but 
a substantially large value of {ip^il'} is favored by light 
quarks existing in the chiral symmetric phase. In this 
sense, practically, one can identify the quarkyonic phase 
as an exotic state where chiral symmetry is restored 
first {{•ipip) ^ 0) and still the confining property remains 
($ ~ 0). In other words the bulk pressure is mostly dom- 
inated by light quarks and, nevertheless, excited quarks 
on top of the Fermi sphere feel a confining force. [There 
is an argument that the confinin g fo rce might cause in- 
homogeneous chiral condensation [1J| . Such a possibility 
is beyond our current scope.] 

There is no strong evidence for such an exotic win- 
dow. Phenomen olog ical considerations could lead to a 
different scenario [15 [^though some suggested arguments 
have been reported 3 and some model studies are 
supportive general the PNJL and PQM mod- 

els rather favor the quarkyonic picture; the model pre- 
dicts the deconfinement temperature weakly dependent 
on fiB- The Polyakov loop tends to be small for any /ie 
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FIG. 1. Chemical freeze-out points taken from Refs. [19l.[20l|. 
The red and blue (upper and lower) bands represent the re- 
gions where the entropy density s and the baryon number 
density n, respectively, increase quickly from 0.3 to 0.8 in the 
unit of free quark-gluon values, Sfrec and nfroc (see Eq. ([T])). 



as long as T is vanishingly small, whereas the chiral con- 
densate melts at high /ub. However, the serious problem 
in any model studies is that the model-parameter choice 
is largely uncertain. The PNJL and PQM models are not 
exceptions. The situation is worse at higher /ib because 
the lattice-QCD data is unavailable then. It should fa- 
tally depends on model assumptions whether the phase 
diagram has the critical point (s) and/or quarkyonic mat- 
ter or even nothing at all. To make any solid statement, it 
is indispensable to impose some constraints on the effec- 
tive model. In this work we attempt to deduce the phase 
structure from the phenomenological point of view. 

Thermodynamics from the Statistical Model Regard- 
ing the QCD phase diagram at finite T and /xb useful 
information is quite limited. Only the chemical freeze- 
out points in the heavy-ion collisions are experimental 
hints about the phase diagram. Although the freeze-out 
points shape an intriguing curve on the fiB-T plane, as 
plotted by error-bar dots in Fig. [U one should carefully 
interpret it. 

The freeze-out points are not the raw experimen- 
tal data but an interpretation through the Statistical 
Model [H, 1131. In view of the fact that the Statistical 
Model is such successful to fit various particle ratios with 
Hb and T only (hq, /is, and /ic are determined by the ini- 
tial condition) , it should be legitimate to take the freeze- 
out points for experimental data, which in turn validates 
the Statistical Model (though why it works lacks for an 
explanation from QCD). 

Let us proceed further accepting that the Statistical 
Model is a valid description of the state of matter un- 
til the freeze-out curve or slightly above. It is then a 
straightforward application of the Statistical Model to 



estimate thermodynamic quantities such as the pressure 
p, the entropy density s, the baryon number density n, 
etc. We shall utilize the open code THERMUS ver.2.1 
to calculate s and n at various T and /ib W^- 

Figure [T] shows the chemical freeze-out points taken 



from Refs. 



on which s and n are overlaid. For 



convenience we normalized these quantities by 
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These are the entropy density and the baryon number 
density of free massless iV^ — 1 gluons and NcNf quarks. 

Here we note that, in drawing Fig.[Tl we have intention- 
ally relaxed the neutrality conditions for electric charge 
and heavy flavors and simply set /ig = fig = t^c — 0. We 
have done so to make it possible to compare the results 
from the Statistical Model to the chiral effective model 
approach in later discussions. [We note that one can 
force the chiral model to satisfy neutrality but it would 
be technically involved [12].] Nevertheless, we would em- 
phasize that the neutrality conditions have only minor 
effects on the bulk thermodynamics and make only small 
differences in any case. 

We should mention that we used Eq. ([T]) with Nc = 
A'f = 3. The choice of strcc and nhcc is arbitrary and 
the following discussions do not rely on this choice, for 
we will use strcc and ntroc just as common denominators 
to display the Statistical Model and the PNJL model 
results. 

The Statistical Model cannot tell us about the QCD 
phase transitions. Still, Fig. [1] is suggestive enough. We 
can clearly see the thermodynamic quantities from the 
Statistical Model blowing up in a relatively narrow re- 
gion. The red and blue (upper and lower) bands indicate 
the regions where s/sf^co and n/uhcc , respectively, grow 
quickly from 0.3 to 0.8. In the Hagedorn's picture 
this rapid and simultaneous rise in s and n has a natu- 
ral interpretation as the Hagedorn limiting temperature, 
above which color degrees of freedom is liberated, i.e. 
color deconfinement. 

Thermodynamics from the PNJL Model Figure [T] is 
useful to have a guess-estimate about the deconfinement 
boundary but we can deduce no information about the 
chiral property. So, to address the QCD phase transi- 
tions, we must find a way to connect the thermodynam- 
ics in Fig. [1] to the order parameters {■iptp) and Here 
let us go into details of the chiral effective model for that 
purpose. 

It is essential to adopt the Polyakov-loop augmented 
model here because the entropy density should contain 
the contribution from gluons which is taken care of by the 
Polyakov loop potential. The PNJL model that we use 
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FIG. 2. Entropy density normalized by SfVee (from 0.3 to 0.8) 
in the Statistical Model (bottom band with red color; same 
as shown in Fig. [TJ and that in the PNJL model with a choice 
To = 200 MeV (top band with green color). The blue band 
between two represents the results with the ansatz ((3}. 



in what follows is defined with the following potential; 

.4 1 
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with a{T) = ao + ai{To/T) + a2{To/T)^ and b{T) = 
b^iTo/T)^. There are thus five parameters one out 
of which is constrained by the Stefan-Boltzmann limit. 
These parameters are fixed by the pure-gluonic lattice 
data as uq = 3.51, ai = —2.47, 02 = 15.2, 63 = —1.75 Q, 
and To = 270 MeV from the deconfinement temperature 
of first order in the pure-gluonic theory. It is important 
to note that only To is a parameter with the mass dimen- 
sion, so that the energy scale is set by this Tq. 

In addition the NJL sector of the PNJL model has five 
more parameters in the three-flavor case [HI; the light 
and heavy quark masses mud and Wg, the momentum 
cutoff A, the four-fermionic interaction strength gs, and 
the U(l)A-breaking six-fermionic interaction strength g^, 
which are all fixed by the pion mass m,r, the kaon mass 
rriK, the eta-prime mass m,,/, the pion decay constant /tt, 
and the chiral condensate ('0^/') [24| . 

In the presence of dynamical quarks, if we keep using 
To = 270 MeV, the simultaneous crossover temperature 
of deconfinement and chiral restoration is over 200 MeV, 
which is too high as compared to the lattice-QCD value. 
It is argued in Ref. that the back-reaction from quark 
loops affects the mass scale To which changes from Tq ~ 
270 MeV for A^f = reduced down to To = 208 MeV for 
N{ = 2 and To = 187 MeV for A^f = 2 + 1 7]. In this 
work we choose to use To = 200 MeV throughout. 

In Fig. [2] we show the entropy density calculated in the 



mean-field PNJL model with To = 200 MeV in the same 
way as presented in Fig. [1] The bottom (top) band in 
red (green) color is the result from the Statistical Model 
(PNJL model). From the figure it is obvious that the 
PNJL model is not consistent with the Statistical Model 
even qualitatively. With the properly scaled To from 
270 MeV down to 200 MeV, the blow-up behavior in s 
from the Statistical Model could be smoothly connected 
to the PNJL model description only for ^ 400 MeV. 
The curvature of the band as a function of /ib is so dif- 
ferent; the PNJL model result is too fiat horizontally and 
it eventually take apart from the region where the Sta- 
tistical Model breaks down. 

Problem and Solution Such a manifest discrepancy 
from the Statistical Model is a crucial drawback of the 
PNJL model. To make the entropy density at /zb > 
500 MeV get saturated, quark degrees of freedom must 
be released at smaller temperature than predicted by the 
PNJL model. 

One can imagine why this happens in the following 
way; the energy scale in the Polyakov loop potential is 
specified by the parameter Tq that may differ depending 
on the effects of T and /ie in the quark sector. We have 
shifted To from 270 MeV down to 200 MeV, through 
which we have incorporated the scale change induced by 
Nf quarks at finite T. In this way we may well consider 
that To should decrease with increasing hb Q- 

Our idea proposed here is to make use of Fig. [5] to 
fix To(/xb) for consistency with phenomenology. [One 
can pick up other thermodynamic quantities than the 
entropy density like the internal energy density, which 
would anyway make little change in the final result.] In 
Ref. [1^ the freeze-out curve is parametrized as Tt(/iB) = 
a — — c/^B with the fitting result a = 166(2) MeV, b = 
1.39(16) X 10-^ MeV"\ and c = 5.3(21) x 10"" McV"^. 
Because the behavior of the entropy density must be 
dominantly controlled by deconfinement, we postulate 
that the change in Tq is to be correlated with Tf(//B)- 
[We see that the freeze-out points and the entropy band 
in Fig. [T] have roughly same curvature.] Let us simply 
use same b and make an ansatz as 

ToiflB 



To 



1 



(3) 



which yields the blue band in the middle of Fig. [2j [To 
prevent unphysical negative Tq for large ^b we set a 
threshold at 10 MeV so that Tq > 10 MeV. Hence, the 
results at T < 10 MeV are not meaningful.] We see at a 
glance that the results from this modified PNJL model 
have a reasonable overlap with the Statistical Model re- 
sults in the whole density region as plotted. 

At this point one might have thought that the energy 
scale in the quark (NJL) sector should be modified as 
well. We will come back to this question after discussing 
the phase diagram next. 
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FIG. 3. Phase boundaries associated with deconfinement 
(blue band) and chiral restoration (red band). Each band 
represents a region where the (normalized) order parameter 
develops from 0.4 to 0.6. 

Phase Diagram Now we get ready to proceed to the 
possible QCD phase diagram that is at least consistent 
with the Statistical Model outputs in Fig. [T] Using 
the standard computational procedure of the mean-field 
PNJL model we can solve {iptp) and $ as functions of 
T and fj,B, from which the phase boundaries of chiral 
restoration and deconfinement are located. 

Figure |3] shows the phase diagram from the modified 
PNJL model. The blue (red) band is a region where the 
Polyakov loop (normalized light-quark chiral condensate) 
increases from 0.4 to 0.6. In contrast to the old PNJL 
model ones, the new results show that the chiral phase 
transition is almost parallel to and entirely above the de- 
confinement, which agrees with the situation considered 
recently in Ref. . We have found the critical point 
at (/iB,T) ~ (45 MeV, 330 MeV), but would not take it 
seriously since its location is easily affected ^5]. Still, it is 
a good news for the critical point search that two QCD 
phase transitions stay close to each other, for the exper- 
imental signature would be detectable only if the critical 
point sits sufficiently near the freeze-out point. 

Discussions It is an intriguing observation that the 
chiral phase transition occurs later than deconfinement. 
This is quite consistent with the Statistical Model as- 
sumption. In the Statistical Model the hadron masses 
are just the vacuum values and any hadron mass/width 
modifications are neglected, which would be a reasonable 
treatment only if the chiral phase transition is above the 
Hagedorn temperature. Under such a phase structure, 
besides, our assumption of neglecting /ie-dependence in 
the NJL-model parameters turns out to be as acceptable 
as the Statistical Model treatment. This can be under- 
stood from the fact that the NJL part yields the hadron 
masses in the vacuum which are intact in the Statistical 



Model. 

The failure of the standard PNJL model is attributed 
to baryonic degrees of freedom missing in a non-confining 
quark description. Hence, one may say that a modifica- 
tion made in stems from such crossover between 
baryons and quarks, which is presumably parametrized 
by the Polyakov loop alone, similarly to the transverse 
gluon pressure. It is an important question how our phe- 
nomenological input (jS]) is validated / invalidated from the 
first-principle QCD calculation, which will be answered 
by future developments in the functional renormalization 
group method [3] . 

Finally, our conclusion is that, if quarkyonic matter is 
defined by restored chiral symmetry with confinement, it 
is not consistent with the Statistical Model and is unlikely 
to occur. However, to complete our analysis it should be 
necessary to think of the quarkyonic spiral |14l] , which is 
an important future problem. 
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